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Abstract. We study the connection between the BCS pairing model and the 
inhomogeneous vertex model. The two spectral problems coincide in the quasi- 
classical limit of the off-shell Bethe Ansatz of the disordered six vertex model. 
The latter problem is transformed into an auxiliary spectral problem which 
corresponds to the diagonalization of the integrals of motion of the BCS model. 
A generating functional whose quasi classical expansion leads to the constants of 
motion of the BCS model and in particular the Hamiltonian, is identified. 
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1. Introduction 

One of the most successful models of interacting electrons is the BCS model of 
pairing Q. Originally proposed to describe properties of superconductors [g], the 
pairing idea has been applied to a large variety of physical systems in nuclear physics || 
and in QCD Recent experiments in metallic nanoparticles |5| have renewed 

the interest in the problem of pairing correlations in mesoscopic systems ||. The 
BCS Hamiltonian consists of a kinetic and an interaction term which describes the 
attraction between electrons in time reversed states 
si si 

H = £ i° c \o c i<y c ]T c ]i c J'J- c J'T ■ (1) 

The quantum number j 6 {1, . . . , f2}, a 6 {|, J,} labels a shell of doubly degenerate 

time reversed states of energy ej ; Cj^ a and cj a are the corresponding electronic 

operators; g is the BCS coupling constant. The low energy properties associated 

to this model are universal functions of the ratio between two energies, the single 

particle average level spacing and the BCS gap Q. 

Various exact results have been obtained for the BCS Hamiltonian. 

In the limit g — > oo the exact eigenvalues and eigenstates can be found (see for example 

Ref. J?]]) and the integrals of motion are Gaudin Hamiltonians ||. An important 

consequence of the relation with the isotropic Gaudin magnet (discussed in Appendix 

A), is that the Quantum Inverse Scattering Method (QISM) |9J for the g — > 00 BCS 
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model can be related with the QISM for the Gaudin model |Tj|. The g — > oo BCS 
model can be also related to the inhomogeneous vertex models~[^2|, |l3) . 
Much less work has been done for finite g. The exact solution was found by Richardson 
and Sherman (RS) |19 and independently by Gaudin by means of the Bethe 
Ansatz (BA) technique. Approximate expression of correlation functions were found 
in the Ref. plf . More recently, the integrals of motion of the BCS model were 
obtained [^2], [ll| and were diagonalized by means of the algebraic BA. 
In this work we show that for finite g the BCS model is connected to a disordered 
six vertex model through the Off-Shell BA (OSBA) introduced by Babujian et al. in 
Refs. J!?], |l^|. In this framework, the known connection between the isotropic Gaudin 
models and the inhomogeneous vertex model is obtained as the Mass-Shell limit which 
corresponds to g — ► oo. 

A strong hint towards our result is provided by a recent work by Sierra |23f| who 
has shown connection between the BCS pairing model and a su(2) c Wess-Zumino- 
Novikov-Witten Conformal Field Theory (CFT), in the singular limit when the 
central charge is infinite; the RS wave functions solve the Knizhnik-Zamolodchikov 
equations for the CFT correlation functions. The results of Sierra are indeed related 
to the connection existing between models in statistical mechanics and correlations 
functions of a suitable CFT established through the OSBA. In fact the solution of the 
quasi classical OSBA equations is equivalent to solve the Knizhnik-Zamolodchikov 
equations |24|, and in particular the quasi-classical OSBA equations for the vertex 
models generate the correlators of the su(2) Wess-Zumino-Novikov-Witten CFT. 

The paper is organized as follows. In section II we review the exact solution of 
the pairing model. In section III the quasi-classical expansion of the OSBA of the 
disordered six vertex model is identified as the diagonalization of the BCS model. 
Section IV is devoted to the conclusions. The connection between the diagonalization 
of the pairing model for infinite pairing coupling constant g and the diagonalization 
of Gaudin magnet is reviewed in Appendix A. In the Appendix B we summarize the 
QISM of the inhomogeneous vertex model. 



2. The exact solution of the pairing Hamiltonian 

In this section we review the exact solution ||, |(J of the BCS model (Eq. (@)) and 
the formulation of its integrability. Due to the form of the pairing interaction in 
Eq (|l|), single occupied states are frozen and we can focus on scattering of pairs. The 
Schrodinger equation for a state of N Cooper pairs 

H\N) = £\N) ; (2) 

has the solution [2(J 



W)= J! ^+(^,^10) ; a + (e a ,e) := ]T • 



l 2£ 3 

N 



(3) 



The operators aj := Cj^Cj^, <r+ = {cjy and cr| := (cJ- ^Cj,^ + c - ^Cj t i — l)/2 realize 
su(2) in the lowest representation. The vacuum state is the highest weight vector of 
su(2): |0) := 1 1/2, — 1/2). The operators {tr ± (e Q , e), a z {e ai e)} generate (for generic 
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e a ) the Gaudin algebra £[sZ(2)] (see the Appendix A and Eqs. (|15|) of the next section). 
The energy £ is given in terms of the spectral parameters e a which satisfy the algebraic 
equations Q 

1 A 2 A 1 

^ " x -=Q, a = l,...,N. (4) 



9 jri ep - e a 2 £j - e 

The method employed by RS has analogies with the coordinate BA technique. In fact, 
in the coordinate BA the ansatz functions are plane waves (describing free particles) 
modified to include the interaction. In the RS solution the ansatz functions are the 
solutions of the model when pairs of time-reversed electrons are treated as bosons; 
these functions are modified because Cooper pairs behave as hard-core bosons. In 
both the RS and the BA procedures the modification enters the set of the algebraic 
equations for the rapidities (Bethe equations) parameterizing the eigenvalues of the 
Hamiltonian. 

By using the spin realization of pair operators {cr^ , c^}, the pairing Hamiltonian 
can be written as a quantum spin model with long range interaction in a non uniform 
fictitious magnetic field, given by £j 

n Q 

H = Y, 2e W - f Yl ( a ? a J + a t a i ) + const ■ ( 5 ) 

3=1 j,l=l 

Cambiaggio et. al. [E2| found that the integrals of motion tj of this model, if Ej ^ £/, 
Vj 7^ I have the form 

Ti = - g a] - 3j , (6) 

and satisfy the commutation relations [H,Tj] — [Tj,Ti] — 0, Vj,l € {l,...,f2}. The 
operators Ej in Eqs. (||) are spin-1/2 Gaudin Hamiltonians § 

1=1 b 3 tl 

The commuting operators tj were also found by Sklyanin [pd]| by taking the quasi- 
classical limit of the monodromy matrix of the inhomogeneous vertex model twisted by 
a term proportional to cr| / g. The pairing Hamiltonian can be expressed as function 
of the integrals of motion as 

- n si 

-3 H = — 2£ j T J + T 3 T i + const - ( 8 ) 

9 9 3=1 3,1=1 

In the limit g — > 00 the problem is equivalent to the diagonalization of (all) the Gaudin 
Hamiltonians (see the Appendix A for details). 

3. The OSBA of the inhomogeneous vertex model and the pairing model 

Vertex models are two dimensional classical models which were solved long ago by 



inverse methods [12 . Generalizations to su(2) higher representations and to include 
disorder were intensively studied |l4|, [l^, [l(| . 

In this section we introduce a vertex model in which the inhomogeneity is due to 
the combination of given (see below) disordered distribution of both the spin and the 
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Figure 1. The vertex (a,j) £ A. Disorder Zj is distributed along the vertical 
lines. It is due to both the spin inhomogeneity and impurities distributed in 
the lattice: Zj = z{Sj) + z(Ij) 



impurities in the lattice. Then we apply the scheme developed by Babujian et. al. [|l7| 
to relate this inhomogeneous vertex models to the BCS model Eq. (|l|). 
The model is defined in the following way. On the edges of the square lattice 
A : N v x Nh (N v columns and Nh rows) are arranged N v + 1 types of spin variables. 
On horizontal edges (labelled by a = 1 . . . Nh) are arranged the spins a taking spin 
projection m a 6 {±1/2}. On the columns (labelled by j = 1 . . . N v ) the spin variables 
Sj can take any value rrij <E {— Sj, . . . , +Sj} of the Sj-th representation of su(2). The 
partition function is restricted to configurations for which an even number of spins 
are into (or out of) each lattice site (vertex); configurations in which the four spins 
are all in or all out are excluded (ice rule). The "scattering" between spin states 
(m a ,mj) — > (m' a ,m'j) of vertex (a, j) (see Fig. have weights fixed by the universal 

matrix elements R™ 3 ' m 3 / (A — zf) where A is the spectral parameter. The quantities 
Zj (see also Ref. fl5|| ) shift spectral parameters as inhomogeneities which we assume 
distributed only along the columns of A (Fig. |l|). 

The variables Zj take into account of disorder induced by the mixture of spin 
representations and/or by the actual distribution of impurities Ij in A: Zj = 
z(Sj ) + z (Ij). We assume that both z(Sj) and z(Ij) enter the universal matrix (see 
Eq. (B.l)) in the same functional form. The disordered six vertex model corresponds 
to the choice: z(Sj) = 0, z(Ij) 0. We impose periodic boundary conditions. 
The transfer matrix T(A|z), where z := [z\ . . . z ^ h ), can be expressed in terms 
of rational i?-matrices Rx, X = {a, S} (see Eq. (B.l)) which fulfill Yang Baxter 
relations (see the Appendix Appendix B| ) . This implies the integrability of the model: 
[r(A|z),T(/i|z)]=0. 

The application of the transfer matrix to the Bethe vector $(Ai . . . Ajy|z) reads: 

T(A|z)$(Ai ...\ N \z)= A(A, Ai . . . X N \z)<!>(X 1 . . . X N \z) 
N p 

-$ Q (Ai..A Q _i, A, A q+ i..Aat|z) , (9) 



, A - A« 



(for the explicit form of the quantities T(A|z), $(Ai . . . Ajv|z), A(A, Ai . . . Ajv|z), 
F a , and $ a (Ai ..A Q _i, A, A a +i..Ajv|z) in @, see Ref. @). The condition for the 
diagonalization of T (and of the constants of motion generated by T) is that 
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the spectral parameters are chosen to cancel the "unwanted terms" (the second 
contribution to Eq. (g)) in the spectral problem (|9|); a sufficient condition is: F a = 
(algebraic BA equations). Such a condition has been termed as "Mass-Shell" 
constraint jlj imposed on Eq. (^|). The OSBA spectral problem, instead, arises when 
the "unwanted terms" are considered in Eq. (||); the spectral parameters obey a new 
set of equations called OSBA equations (see Eq. ( fL3| ) below), unfixed). The statement 
of the OSBA problem has been shown remarkable in the quasi-classical limit. In fact, 
it was discovered that the solutions of the quasi-classical OSBA problem satisfy the 
Knizhnik-Zamolodchikov equations J24|, |25| for the su(2) CFT. In the following we 
shall see how the quasi-classical limit of the OSBA problem for the disordered vertex 
model is solved by spectral parameters fulfilling the Eqs. (jj). 

The quasi-classical limit of the vertex model is obtained through the expansion 
of Rx(x; rf) in powers of 77 (Rx{x; 0) is the identity). 



Using the expressions for the monodromy, transfer, and universal matrices Eqs. (B.l), 



(B.3), the quasi-classical limit of the OSBA Eq. (|9j) reads 

j—l a—l 

up to order 0(rj N+2 ), where the explicit form of h and (j> a in ( JlCj ) is given in the 
Ref. 0. By integrating Eq. (|l(]) on a closed loop in the complex A-plane encircling 
the pole A = Zj we obtain 

N 

^ = m+Ettt s M' (u) 



— Zi 
= 1 J 



a=l 



N v N 
SlSj 



1 = 1 * a=l 
N N v 

/« = Er- L T-Eir-' a = l...N. (13) 
The Bethe vectors in the quasi-classical limit are 

N N 

:= [] S+(X , z)\vac) , 4>' a := JJ S+(A/3, z)|mc) , (14) 

Here |uac) = <8>j=i|sj, — Sj), where S^~|sj, — Sj) = 0, i.e. \vac) is the highest weight 

vector in ®jsu(2)j. The three operators S ±,z (Xp, z) := YljZi Sf' z /(zj — A/3) generate 
higher dimensional representations of the Gaudin algebra Q[sl{2)], given by (see also 



Eq. (A.2)) 



[s*(x a ,z),s±(x 0> z)] =± s ± {* a ,*)-sr k (>»,*) 

Ap — A a 

[S+(X a , z), S-(X , z)] = S'M-S'(X p ,z) 

A/3 — A Q 

The "Mass-Shell" constraint: f a = corresponds to the diagonalization of the 
Gaudin model (see Appendix A). 
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The solution of the spectral problem for the pairing model is recovered 
substituting 



a = l...N, (15) 




in the left hand side of Eqs. (|l3| ). In fact, the resulting equations coincide with Eqs. (|j). 
Substituting Eqs. ( |l5| ) in (|ll[), and summing over index j = 1 . . . N v we obtain 

= 0, (16) 

where we have used the fact that: Ylj=i hj — ~ Ea=i Sj=i s j/{^a ~ Zj). Eq. ( [l6| ) 
shows that, the OSBA spectral problem is transformed in a spectral problem involving 
only diagonal matrix elements of suitably shifted (by f a ) transfer matrix of the vertex 
model (in the quasi classical limit). 

Since the limit g — > oo should correspond to the same result of fa ^ for generic Sj 
(compare with Eqs (A.4)), we impose that the distribution of spins Sj through the 
lattice fulfills the condition 

N ", n -9.*a 

. (17) 



(1 - 2S j) _ 

p[ (^-Zj) 



In this case Eq. ( |l5| ) reduces to 

fa = ^~, a=l...N. (18) 

We choose Sj = 1/2 Vj in order to fulfill Eq. @: the inhomogeneous vertex model 
becomes the disordered six vertex model since z(Sj) = and z(lA 7^ 0. This implies 
that 

Hj = Sj , and <j> = \N) , 

where: N v = O (compare with Eqs @, (@)). Eq. @ can be transformed in the 
following eigenvalue equation: 

UN/ _ (jZ y \ n N 



1 1 , 



n 9 N">J> (20) 

where Eqs. (§ (or ©, © ), ©, and: E?=i **(e«, e)^ - 1/2 E?=i E?=i V(e« - 
2>£j)4> have been used (the parameters in Eqs. (p~2[ ) — (|20|) are redefined as Zj <-> 
2s j , A Q <-> e Q ) . We point out that quantities in fl20|) are the eigenvalues of operators 
7j in Eq. (g) for generic fi/JV. At "half filling" ft = 2iV Eq. @ reduces to 

-('- 

Equations (pcf), ( pl| ) coincide with those ones found by Sklyanin |Ll| and by Sierra p3| |. 
The main result obtained in this paper is the connection between Eqs. (^Tj) and 
Eqs. (19). (p0|) through the Eq. (fl5|). The OSBA problem for the disordered six vertex 
model (which does not account for diagonalizing the transfer matrix of the vertex 



,.. 77 " 2/ , (j = l,. ..,«)• (21) 
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model) reveals the existence of a class of spectral problems (parameterized by f a ) 
which turns out to be diagonal on the quasi-classical Bethe vectors basis. For f a fixed 
by Eqs. ([l8]) the diagonalization of the BCS model is obtained. 

Furthermore, what we have discussed so far implies that the pairing Hamiltonian can 
be obtained from functionals of Tj whose quasi-classical expansions have the following 
form: 



T(e|z)=^r ? 2a e°- 1 

a=0 



with 



r(e) 



12 

E 



2e,- 



(22) 



(23) 



We point out that [T(e|z), T(e'\z)] = , Ve, e' since quantities Tj's commute each 
other. The residue in the poles e = 2ej of the rj 2 coefficient provides the integrals of 
motion Tj. The residue of the rf coefficient reads (see Eq. (||)) 



n 

E 

3,1=1 



Tin 



l - l 



(24) 



4. Conclusions 



We have established a novel connection between the disordered six vertex model 
and the BCS model for generic g through the OSBA procedure. The BCS model 
is diagonalized by the quasi classical limit of the OSBA equations of the disordered 
six vertex model. Retaining certain off diagonal terms of the transfer matrix of the 
vertex model corresponds to the diagonalization of the integrals of motion of the 
pairing model for finite g. The "mass shell" condition (and then the diagonalization 
of the quasi classical transfer matrix of the vertex model) reproduces the limit g — > oo; 
the corresponding problem is the Gaudin spectral problem. 

The integrals of motion of the BCS model coincide with the integrals found by 
SklyaninJTl| by considering a twist in the monodromy matrix of the vertex model 
(see Eq. |B.3| ). The algebraic equations which diagonalize these integrals of motion via 
algebraic BA (namely via the Mass Shell BA procedure) coincide with the Richardson's 
equations. The present study shows that the Sklyanin's procedure produces the same 
results of the Off Shell BA procedure applied to the untwisted monodromy matrix. 

The existence of the relation between BCS model and quasi-classical vertex 
models, found in the present paper, is consistent with the correspondence between 
CFT and the BCS model recently found by Sierra 

Eqs. (||) were already conjectured by Gaudin (see formulas (5.15), (5.16) of 
Ref. ||) as connected (through Jacobian of certain matrices) with the norms of the 
Bethe vectors, det(df a /dep) ~ ||<?!>|| (see Refs. plj , p6fl). In this work we have shown 
that the Jacobian is connected with OSBA of the vertex model. This might be useful 
to compute norms (and scalar products) and, then, to express the correlation functions 
of the BCS model as suitable determinants. This exact calculation is our major task 
in the future. 
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Appendix A. The pairing model and the Gaudin spectral problem 

In this Appendix we discuss the connection between the pairing model and the Gaudin 
model. 

The limit g — » oo of the constants of motion tj (^) coincides with Hamiltonians Sj . 
Since Eq. (||), the spectrum of the pairing problem coincide with that one of the 
Gaudin magnet [||: S(u) := Ylj=i^j /( u ~ ^ £ j) ( u * s a complex parameter). The 
total energy is: h(u) :— Y^j=i hj/(u — 2ef) (hj is fixed by Eq. ([l2|)). The Bethe 
vectors of the Gaudin and the pairing problems coincide formally for any g since 
operators (cr ± (u, e), cr z (u, e)) in (|^) generate the Gaudin algebra C/[s^(2)] in the lowest 
representation: 

[a z {u,e),a ± (w,e)} = ± y — L - L , (A.l) 



[cr + (w,e),er (w,e)} 



w — u 
a z (u, e) — <t z (w, e) 



where a z (u, e) := Y,j=i a j/ ( 2e j _ u )- 
However, the spectral parameters entering the eigenvectors of the two models satisfy, 
for generic g, different equations (compare with Eqs. 
n 2 n 1 



We point out that Eqs. ( |A.2| ) are the limit g —> oo of Eqs. (Q) for the pairing model. 
In this limit the two models have the same eigenvectors (see Eqs (||), (0), (A. 2)). 



Thus the diagonalization of the Gaudin model is equivalent to the diagonalization of 
the BCS model for g — > oo. 
The limit of large g in Eq. (|§|) gives: 



■ i i 3 tl 



£ 3< J j G l 



= -f E a t °i < ( A - 3 ) 

J,/=l 

which (consistently) reproduce the Hamiltonian ([!]) for large <?. 

The QISM was applied to Gaudin magnet for generic spin Sj in Ref. ^ |27|| . In 
this case the Gaudin Hamiltonians are Hi :— ^? =1 • Sj/(ej — £/). The spectral 
parameters obey 
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where — Sj corresponds to the highest weight vector of Sj. 

Appendix B. Integrability of the inhomogeneous vertex model 



In this appendix we discuss the QISM of the inhomogeneous vertex model together 
with its quasi-classical expansion. 

The universal matrix of the model reads 

R X (X - z; r?) = 1 ® 1 + /(A - z, f])(J ® X, (B.l) 

where f(x,rf) :— 2r\j(r] — 2x) depending on the arbitrary parameter 77 6 M. The 
physical i?-matrix of the model corresponds to take X = S in (B.l); the auxiliary one 
corresponds to X = a and z = 0. Both these matrices fulfill homogeneous (S = a) 
Yang-Baxter relations; in addition they fulfill the the inhomogeneous (S ^ a) one Q|: 



R?(\-ti)Rg(\-z)Rg(ji-z) = 

Rf(»-z)Rf(\-z)R?(\- 

where R 12 = R® 1 and i? 23 = 1®R act on the vector space V\ <g> V2 



(B.2) 

3V3; R 12 and R 23 



act as the identity in the vector spaces 3 and 1 respectively. Eq. (B.2) is the sufficient 
condition for which the model, can be solved by diagonalizing the column to column 
(which is 2 x 2 matrix in operators Sj) transfer matrix (instead of the row to row 
one which is (2s j + 1) x (2s j + 1)) obtained through the trace in the two dimensional 
horizontal vector space (which is spanned by spins along the rows of A) labelled by 
"(0)": 



T(X\z) :=tr i0) J(X\z) 



(B.3) 



Twisted monodromy matrix is obtained by J(A|z) — > e OCT J ; J(A|z). The transfer 
matrices commute at different values of spectral parameters: [T(A|z), T(fi\z)] = 
(it proves the integrability of the model) since the monodromy matrix J(A|z) := 



■ zj) satisfies: R (T (X — fi)[J(\\z) ® J(/x|z)] = {J(fi\z)® J(X\z)]R a (X — (i) 



II] 

(induced by Eq. (B.2). The matrices R^ and R a arc 

l + f(x,7 1 )S 
f(x,v)S+ 



R^(x, V ) := 



f(x,r))Sj 

l-f(x,Tj)S* 



R a (X - pb,rj) := 



( 1 



V 







1 / 



(B.4) 



(B.5) 



where: b(X — jj,) := 77/(77 — A — /i), c(X — fi) := X — /j/(A — /i — 77) (note that R r7 (X — /1, 77) 
is z-independent). 

In the quasi-classical limit, the system generates a hierarchy of integrable systems 
in the quasi-classical limit since 



E 1 

a=b+c=0 



T b (X\z),T c (fi\z)}=0, 



(B.6) 
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where we have used the ^-expansion of the transfer matrix: T(A|z) = X^^Ln V a Ta,(M z ) 



(the sum in Eq. (B.6) is meant on ordered partitions of a including b V c = 0). Up to 
order ij 2 , the transfer matrix reads 

N v 

T(A|z)=21 + 2 ?7 2 ^— (B.7) 

j=i Z ' } 



where the Hamiltonians Hj in Eq. (B.7) are 



N v 



Hj = j2—> (i = i,...,i\y. (b.. 

L ^ y. — 



1=1 

1*3 



Zl 



The transfer matrix (B.7) coincide with the quasi-classical expansion of the twisted 



Gaudin model's transfer matrix (see formula (1.16) of the Rcf. [|ll| 
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